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ABSTRACT 

The transient X-ray source CX0M31 004205.77+411330.43 exhibited several out¬ 
bursts during our long-term monitoring campaign of ^monthly observations of the 
M31 center with Chandra. However, the decay profile appears to be unlike those ob¬ 
served from Galactic transients. We followed up the 2011 outburst with two band 
HST /ACS observations, one in outburst and the other in quiescence, and used differ¬ 
ence imaging to search for a counterpart; this would be dominated by re-processed 
X-ray emission from the disk. We found a counterpart with B = 28.21+0.16. An XMM- 
Newton observation from a previous outburst yielded a spectrum that is well described 
by an absorbed power law with absorption equivalent to ~2.6xl0 21 H atom cm' 2 and 
photon index ~1.8; the highest quality Chandra spectrum yielded only ~130 counts, 
and best fits consistent with the XMM-Newton results. We calculated an absolute V 
magnitude of +1.9 during outburst for a typical disk spectrum. An empirical relation 
between the ratio of X-ray to optical flux and orbital period suggests a period <4 hr 
for a black hole accretor. Such a short period is expected to yield an asymmetric, pre¬ 
cessing disk, and we propose that the observed decay lightcurve is due to modulation 
of the mass transfer rate due to the disk precession. 

Key words: x-rays: general — x-rays: binaries 


1 INTRODUCTION 

We have been monitoring the central region of M31 with 
Chandra, averaging ~1 observation per month over 1999- 
2012, in order to discover X-ray transients. Since then, we 
have reduced our observation rate to 5 per year. Particu¬ 
larly bright or otherwise interesting transients are followed 
up with two HST ACS observations, the first is taken a 
few weeks after outburst, and the second observation is nor¬ 
mally taken ~6 months later; this allows us to identify the 
coun terpart via difference imaging (see e.g. 1 Barnard et al.1 
l20ll and references within). We summarized the results 
of the first 12 transients (labe led T1-T12) found via this 
effort in iBarnard et al.l d2012l ). and performed this anal¬ 
ysis for two further transients: XMM J004243.6+412519 
(known as M31 ULX2, iBarnard et alj|2013al l. an d CXQM3 1 
004252.457+411631.17 (referred to as T13, IBarnard et al.l 
l2014bl) . 

In this work we report our findings for CXOM31 
004205.77+411330.43, hereafter referred to as T14. We note 
that T14 is source number 71 in our variability survey 
of 528 X-ray sources within ~20 / of the galaxy nucleus 


where we made ^200 new X-ray binary (XB) identifica¬ 
tions (Bar nard et al.lf2014al ): the 0.3-10 keV luminosities 
for these sources ranged over ~l-6000xl0 35 erg s -1 . The 
initial outburst of T14 was in 2011 July; however, the sec¬ 
ond observations was made almost 3 years later, in 2014 
June. 

X-ray transients within the central region of M31 are 
most likely to contain a low mass seconda ry, as t he major¬ 
ity of stars there are old (see e.g. [Williams et al.|[200a ). 
Low mass X-ray binaries (LMXBs) may be transient X-ray 
sources due to instabilities in their accretion disks; the disk 
has two stable phases (hot and cold), and an unstable inter¬ 
mediate phase— matter accumulates in the disk in the cold 
phase, and is rapidly dumped onto the compact object in 
the hot phase (see e.g. lLasotall2001l ). However, the X-rays 
produced by accretion from the hot disk prevent the disk 
from cooling; the X-ray luminosity decays exponentially if 
the whole di sk is ionised, and linearly if only part of the disk 
is ionised rtKing fe Ritteijll998l ). 

Ivan Paradiis fe McClintockl (119941) found an empirical 
relation between the ratio of X-ray and optical luminosities 
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of Galactic X-ray binaries and their orbital periods, sugges¬ 
tive that the optical emission is dominated by reprocessed 
X-rays in the disk; this relation holds over a 10 magnitude 
range in optical luminosity, and appears to be independent 
of inclination. Their chosen X-ray band was 2-10 keV. For 
an irradiated accretion disk with radius a , X-ray luminos¬ 
ity Lx, optical luminosity L opt , and temperature T, T 4 oc 
Lx/a * 2 , while the surface bri ghtness of the disk , S, oc T 2 for 
typical XBs (Ivan Paradiis fc McClintockl 1 19941 ). Since L op t 
oc Sa 2 , L opt oc L^ 2 a; also a oc P 2 /^ 3 , where P or b is the orbital 
period. 

Ivan Paradiis fc McClintockl d 19941 ') defined X = 
(L x /Lbdd ) 1/2 (Porb/lhr) 2/3 , choosing Ledd = 2.5xl0 38 
erg s' 1 as a normalizing constant, and found 

My = 1.57(±0.24) - 2.27(±0.32) log X. (1) 

However. Ivan Paradiis fc McClintockl Jl994h sampled a mix¬ 
ture of neutron star and black hole binaries, in various spec¬ 
tral states. A cleaner sample was obtained by A. Moss et al. 
(2015, in prep), who used only black hole transients at the 
peaks of their outburst, and found 

My = 0.84(±0.30) - 2.36(±0.30) log X. (2) 

We note that these two relations only differ significantly 
in normalization, caused by black hole X-ray binaries hav¬ 
ing larger disks than neutron star binaries with the same 
period. We have period estimates for 14 M31 transients 
(T1-T12, ULX2, and T13) observed by Chandra and HST 
dBarnard et al J[20 1 2l . |20 1 3bl . 120 14 IjI ). These period estimates 
appear to be slightly shorter over all than for Galactic tran- 
sients, ranging from <1 hr to 140^ 3 3 hr, and mostly <30 hr 
dBarnard et al.ll2012l. l2013bl. l2014bl) . 


2 OBSERVATIONS AND DATA REDUCTION 

We observed the region containing T14 in 101 Chandra ob¬ 
servations (60 ACIS, 41 HRC). We observed T14 in out¬ 
burst for 5.2 ks with the HST/ACS on 2011 July 21 us¬ 
ing the F435W filter and the WFC mode (or “aperture”) 
(HI, jbml05010, PI M. Garcia); the 4.8 ks observation in 
quiescence occurred on 2014 June 26 (H2, jc6b02010, PI R. 
Barnard), using the WFC1-CTE aperture. The WFC1-CTE 
aperture improves the charge transfer efficiency by moving 
the target to the corner of a detector, allowing point sources 
to be more circular, but does not effect sensitivity; hence 
the level of background observed in H2 is unaffected by this 
change. We also analyzed a 2001 June 29 ~60 ks XMM- 
Newton observation containing a previous outburst that had 
sufficient counts for spectral fitting (0109270101, PI K. Ma¬ 
son). 

All optical analysis was performed with PC-IRAF Revi¬ 
sion 2.14.1, except where noted. The Chandra and XMM- 
Newton observations were analyzed with CIAO version 4.6.3 
and SAS version 13.0 respectively; X-ray spectra from both 
telescopes were treated with xspec version 12.8.2b. 

2.1 HST analysis 

Each HST observation included four flat-fielded (FLT) im¬ 
ages, and one drizzled (DRZ) image. The flat-fielded images 
are corrected for instrumental effects, but not background 


subtracted; the total number of counts in each pixel is given. 
The nati ve ACS resolution is com parable to the FWHM of 
the PSF dFruchter fc et aPl2009lh The drizzled image com¬ 
bines the flat-fielded images, removes any cosmic rays, and 
subtracts the sky background; it is normalized to give the 
number of counts per second per pixel. We used the DRZ 
images from HI and H2 to create a difference image; how¬ 
ever, we used the H2 FLT images for our aperture photome¬ 
try because the DRZ images sometimes contain pixels with 
slightly negative values, and this can cause problems when 
estimating the number of photons in a region. 

2.1.1 Creating a difference image 

We reprojected the H2 DRZ image into the coordinates of 
the HI DRZ image, to produce an accurate difference image. 
To do this, we first registered the HI and H2 images to the 
LGS Field 5 image with ccmap, using unsaturated stars 
that were close to the target. Then, we used the iraf task 
wregister to make the pixel orientation of the H2 image 
match that of the HI image. We registered the H2 image 
to the HI image before mapping to reduce the noise during 
image subtraction. The difference image was produced by 
subtracting H2 from HI using the ftools task farith. 

2.1.2 Measuring the optical counterpart 

For HI and H2 we found the number of photons for each 
FLT image within various extraction regions, yielding Ctot 
photons in total over T seconds. We subtracted the H2 total 
from the HI total to get Cnet source photons. We converted 
this to Vega B magnitude via 

B ~ -2.5 log [Cnet/T] + ZP, (3) 

having obtained the zero point (ZP = 25.7 7 ) from th e ACS 
Zero Point calculatoiQ: we see from ISirianni et al.l (120051 . 
Equation 12 and Table 18) that the conversion from F435W 
counts to B magnitude is within 3<r of our Z P for B-V = 
—0.09 (assuming a typical disk spectrum iLiu et ahl 1200 li b 
We also estimated the 4cr upper threshold for source detec¬ 
tion from the total photon count from H2 (Ctot.m): 

B ~ -2.5 log [4 x C t ° ot 5 H2 /T] + ZP. (4) 

We can convert from B magnitude to My via 

My = B + 0.09 - A h x (1 + 1/3) /1.8 x 10 21 - 24.47, (5) 

where Ah is the line of sight absorption; this accounts for 
the difference in B and V magnitudes of a typical accretion 
disc, a relationship between B band extinction and measured 
line-of-sight absor ptio n towards t he object, and the distance 
to M31 fsee lBarnard et ad1l2012l . and references within). 

2.2 XMM-Newton analysis 

The XMM-Newton observation suffered substantial flares 
in the particle background due to soft Solar protons. We 
selected good time intervals by making a 10-12 keV pn 

1 http://www.stsci.edu/hst/acs/analysis/zeropoints/zpt.py 
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lightcurve from the whole image with 100 s binning, us¬ 
ing standard filters. We extracted source source and back¬ 
ground spectra from the pn, MOS1 and MOS2 data: source 
spectra from a circular region that was optimized by the soft¬ 
ware (17” radius), and background spectra from a nearby 
circular region on the same CCD with no point source (35” 
radius). We grouped the source spectra to give a minimum 
of 20 counts per bin, created the corresponding response 
matrixes and ancillary response files, then fitted the spectra 
with xspec. The pn, MOS1, and MOS2 spectra were fitted 
simultaneously with constants of normalization to account 
for differences in detector responses. 

We obtained uncertainties in the best fit parameters by 
simulating 1000 spectra with the multifake command, with 
deviations in simulated spectra being drawn from the prop¬ 
erties of the observed spectrum. The best fit parameters for 
these spectra were each ranked in ascending order, with 1 <t 
uncertainties derived from the 160 th and 840 th data points. 

2.3 Chandra Analysis 

2.3.1 Locating the X-ray source 

We used 27 X-ray bright globular clusters (GCs) to regis¬ 
ter a combined ~350 ks ACIS image (supplied by Z. Li) to 
the B band Field 5 image of M 31 provide d by the Local 
Galaxy group Survey (LGS) (jMassev et al.l !2006b We used 
PC-IRAF v2.14.1 to perform the registration, f o llowin g the 
same procedure as described in lBarnard et all (|2012|'I . The 
X-ray and LGS positions of the 27 X-ray bright GCs were 
determined using imcentroid; the equivalent FK5 coordi¬ 
nates were calculated for the X-ray and optical position of 
each GC using xy2sky v2.0, distributed with ftools. The 
X-ray positions of each GC were altered to match the LGS 
positions, allowing the registration of the merged Chandra 
image to Field 5 using the iraf task CCMAP. This registra¬ 
tion yielded la positio n uncertain ties of 0.11” in R.A., and 
0.09” in Dec fearnard et al]l2012h . 

The final uncertainties in the X-ray position of T14 
combine the position uncertainties in the X-ray image, and 
the uncertainties in registering the merged Chandra image 
to the M31 Field 5 LGS image. 

2.3.2 Spectral fitting 

Even the best Chandra ACIS spectrum for T14 yielded few 
counts. We modeled the spectrum in XSPEC in two ways: 
i) using a small number of bins with ^20 counts per bin 
and yf statistics; and ii) using a larger number of bins with 
^5 counts per bin and Cash statistics. We estimated the 
uncertainties for each spectral fit using the same method as 
for our analysis of the XMM-Newton spectrum. 

2.3.3 Estimating luminosities 

For each of our Chandra ACIS observations we extracted 
source and background 0.3-7.0 keV spectra from circular 
regions with 10” radii. We then created a response matrix 
using MKACISRMF, and obtained an ancillary response file 
from mkarf. Spectral analysis was performed using XSPEC, 
giving unabsorbed 0.3-10 keV luminosities for each observa¬ 
tion, allowing us to create a long-term luminosity lightcurve. 


Most of the Chandra observations of T14 yielded too 
few photons for spectral fitting. However, we are able to 
estimate the 0.3-10 keV flux (and therefore luminosity) by 
assuming a particular emission model. If this model is a good 
approximation, then the luminosities of closely-spaced ACIS 
and HRC observations should agree; however, if the model is 
inappropriate, then the ACIS and HRC luminosities should 
be systematically offset, due to differences in instrumental 
response. 

For each faint ACIS spectrum we used XSPEC to esti¬ 
mate the unabsorbed 0.3-10 keV flux equivalent to 1 count 
s -1 , assuming the best fit emission model for the XMM- 
Newton pn spectrum of T14. Multiplying this number by 
the background-subtracted intensity yields a flux that is cor¬ 
rected for instrumental effects, and we o btain the luminos¬ 
ity by assuming a distance of 780 kpc (IStanek fe Garnavichl 
I1998h . It was important to do this for every observation be¬ 
cause the roll angle was not constrained, meaning that the 
source appeared in different parts of the detectors, at various 
off-axis angles. 

For HRC observations, we included only PI channels 
48-293, thereby reducing the instrumental background. We 
used the WebPIMMS tool to find the unabsorbed luminosity 
equivalent to 1 count s _1 , assuming the same emission model 
as for the ACIS observations. We created a 1 keV exposure 
map for each observation, and compared the exposure within 
the source region with that of an identical on-axis region, 
in order to estimate the exposure correction. We multiplied 
the background subtracted, corrected source intensity by the 
correction factor to get the 0.3-10 keV luminosity. 


3 RESULTS 

3.1 Locating and measuring the optical 
counterpart 

The uncertainty in X-ray position for T14 is 0.3” in R.A. 
and 0.4” in Dec llBarnard et alJl2014al). Registering HI to 
the LGS M31 Field 5 image provided bv lMassev et al. (l2006h 
yielded r.m.s. offsets of 0.04” in R.A. and 0.02” in Dec. The 
uncertainties arising from HST are negligible in comparison 
with the X-ray uncertainties. In Figure [T] we present details 
of the HST/ACS WFC F435W images of T14 from HI (left), 
and H2 (middle), along with the difference image (right). 
For each panel, an ellipse represents the 3<r uncertainty in 
X-ray position with respect to the LGS image (0.9” in R.A., 
and 1.2” in Dec.). For the difference image, white stars were 
brighter in HI, while dark stars were brighter in H2. We 
find evidence for a faint optical counterpart at 00:42:05.818 
+41:13:30.00, with respect to the LGS Field 5 B band image. 
This is ~8 ; offset from the ce ntres of the Chandra images. 

In lBarnard et all d2014bi) we used a circle with 3 pixel 
radius (3xFWHM) to estimate the flux for T13; however, 
this was impossible for T14 because this region was severely 
contaminated by cosmic rays in the H2 image, and the H2 in¬ 
tensity (7.6±0.08 count s _1 ) was substantially higher than 
the HI intensity (5.5 count±0.04 s _1 ). Extracting counts 
from a 3x3 pixel region yielded very similar result to ex¬ 
traction from a 2x2 pixel region: background-subtracted in¬ 
tensities of 0.11+0.02 and 0.10+0.02 count s _1 respectively. 
All uncertainties in this work are quoted at the lcr level. 
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Figure 1. Left: Detail of the 2011 July 21 HST/ACS WFC (HI) image of T14 during outburst. Middle: 2014 June (H2) image of T14 in 
quiescence. Right: Difference image, where white points were brighter in HI, and dark points were brighter in H2. North is up, East is 
left. The 3c uncertainty in X-ray source position (0.9” RA, 1.2” Dec) is indicated by an ellipse. We highlight the counterpart with two 
lines. 


For HI, the 3x3 pixel region yielded 6688 photons over 
5240 seconds, while the 2x2 region yielded 3406 photons; 
for H2 the 3x3 region produced 4263 photons over 4790 
seconds, while the 2x2 region gave 2658 photons. These re¬ 
sults allow us to estimate the B magnitude of T14 to be 
28.21±0.16 in HI (using Equation 3) , with a 4er detection 
limit of 28.6 magnitudes measured for H2 from the 3x3 pixel 
region (Equation 4). 

The faintest counterpart that we detected previously 
had an apparent B magnitude of 24.87±0.09, but we mea¬ 
sured 4a upper limits of B >28.7 for two other transients; the 
detection limit is strongly location-dependen t, ranging from 
B ^26 in th e bulge to B ~29 further out llBarnard et al.l 
120121 . l2013bl ). We note that the HI DRZ file, which is 
background-subtracted, yields an intensity of 0.16 count s -1 
for the 3x3 region, and 0.10 count s^ 1 for the 2x2 region, 
which are consistent with our intensities obtained by sub¬ 
tracting H2 from HI. 

3.2 Characterising the X-ray spectrum 

3.2.1 XMM-Newton 

The XMM-Newton pn spectrum of T14 contained ~550 net 
source counts, with an intensity of 0.0209±0.0009 count 
s _1 , while the MOS1 and MOS2 spectra contained 242 and 
323 counts respectively; this is sufficient for simple emis¬ 
sion models. For the simultaneous spectral fits, k\ and 
fc 2 represent the normalizations for MOS1 and MOS2 re¬ 
spectively. X-ray binaries exhibit characteristic emission 
spectra that depend on the accretor and spectral state. 
The hard state is common to all X-ray bina ries, whether 
the a ccretor is a neutron star or black hole (Ivan der Klisl 
1 19941) : its emission is dominated by inverse Comptoniza- 
tion of cool photons on hot electrons, and may be repre- 
sen ted by a power law with photon index (T) = 1.4-2 .1 (see 
e.g. Ivan der Klisll 19941 : IRemillard fc McClintockll2006l ). This 
hard state is obs erved at luminosities below ^10% of th e 
Eddington limit dGladstone et al.l 120071 : iTang et alJ l201ll) . 
Black hole transients often exhibit a thermally dominated 
state during outburst, characterized by a multi-temperature 
disk blac kbo dy with inner disk temperature (feTin) ~1 keV 
(IRemillard fc McClintockil2006l ). The other states of NS and 
BH XBs are more complex, with multiple emission compo¬ 
nents, so we do not attempt to fit them here. 


data and folded model 

pn, MOS1, MOS2 



Energy (keV) 

Figure 2. Unfolded XMM-Newton pn, MOS1 and MOS2 spec¬ 
tra for T14 from the 2001 June 29 observation of T14, fitted 
with the best fit absorbed power law model. The y axis shows 
the model multiplied by channel energy. For this model, JVh = 
2.6±0.3xl0 21 atom cm" 2 and T = 1.84T0.12, x 2 /dof = 42/54; 
the MOS1 and MOS2 spectra were normalized by factors of 
1.02±0.03 and 1.15±0.03 respectively in order to account for dif¬ 
ferences in detector response. 

Fitting an absorbed power law model yielded a best fit 
line of sight absorption (IVh) of 2.6±0.3xl0 21 atom cm~ 2 
and T = 1.84±0.12, fci = 1.03±0.03, fc 2 = 1.15T0.03; the 
0.3-10 keV unabsorbed luminosity was 1.89±0.12x 10 37 erg 
s _1 , and y 2 /dof = 42/54. 

The best fit disk blackbody model yielded IVh = 
0.09±0.03x 10 21 atom cm -2 and fcTin = 1.460±0.13 keV, ki 
= 1.00±0.03, fc 2 = 1.15±0.04; the 0.3-10 keV unabsorbed lu¬ 
minosity was 1.25±0.08 xlO 37 erg s” 1 and y 2 /dof = 56/54. 

The unabsorbed luminosities obtained from these emis¬ 
sion models are inconsistent at the 3 a level; this is due to 
the lower column density from the disk blackbody model, 
caused by the natural downward curvature of the model at 
low energies. The fit may be improved by using a power law 
model rather than a disk blackbody model: Ay; 2 = 7 for 3 
free parameters, and the probability that this improvement 
is genuine is 93%, i.e. not 3<r. When we consider that the lu¬ 
minosity is consistent with the hard state of any XB (NS or 
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Figure 3. Unabsorbed 0.3-10 keV lightcurve for T14, from ACIS 
(filled circles), HRC (hollow circles) and XMM-Newton (trian¬ 
gle) observations, assuming an absorbed power law with TVjj = 
2.6xl0 21 atom cm -2 and V = 1.84. We also show indicate times 
of the HST observations with downward arrows. 

BH accretor), then we find the power law emission model to 
be most likely. We present the unfolded pn spectrum for T14 
assuming the best fit power law emission model in Fig. [2] 
the y-axis shows the spectrum multiplied by channel energy 
in order to show the distribution of flux. 

3.2.2 Chandra 

The best Chandra ACIS observation of T14 yielded ~130 
counts in the 0.3-7.0 keV range, with no photons detected 
in the background region over the ~4 ks exposure time. This 
produced a 6 bin spectrum for y 2 fitting, and a 25 bin spec¬ 
trum for fitting with Cash statistics. 

The best y 2 fit absorbed power law model yielded Ah 
= 8±6xl0 21 atom cm -2 , T = 1.8±0.5, and a 0.3-10 keV 
luminosity of 9±3xl0 37 erg s _1 . The best Cash statistic 
fit yielded Ah = 4±3xl0 21 atom cm -2 , T = 1.5+0.3, and 
0.3-10 keV luminosity 7.1±l.lxl0 37 erg s _1 . These fits are 
consistent with those obtained from the XMM-Newton spec¬ 
trum. 

Fitting an absorbed disk blackbody model with x 2 
statistics resulted in Ah = 3±3xl0 21 atom cm -2 , kT = 
2.1±0.7, and L = 5.2±0.8xl0 37 erg s _1 . Cash statistics 
yielded Ah — 2±2xl0 21 atom cm -2 , and kT = 2.0±0.7, 
with a 0.3-10 keV luminosity of 3.1±0.6xl0 3 ' erg s _1 . 

3.3 The X-ray lightcurve 

We obtained a long-term 0.3-10 keV unabsorbed luminosity 
lightcurve for T14 from the ACIS and HRC observations 
by assuming an absorbed power law emission model with 
Ah = 2.6xl0 21 atom cm -2 and T = 1.84. In Fig. [3] we 
present luminosities from the ACIS (filled circles) and HRC 
(hollow circles) observations; we indicate the times of the 
HST observations with downward arrows, and also include 
the XMM-Newton observation (triangle). 

We see that T14 exhibited several outbursts; they do 
not appear to follow the traditio nal exponential or l inear 
decay seen in Galactic transients (living fc Ritterill99a) : we 
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4500 5000 5500 
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Figure 4. Close-up of the outburst examined by Chandra and 
HST. The axes and symbols are the same as for Fig. [3] 

discuss a possible reason for this below. The luminosities ob¬ 
tained from XMM-Newton, ACIS and HRC are all consistent 
for the first outburst despite the very different instrumental 
responses, and the ACIS and HRC luminosities appear to 
be in good agreement for the whole lightcurve, showing that 
the emission model is probably a good approximation to the 
true spectrum. 

However, we caution that we did not observe any pe¬ 
riod where the luminosity exceeded 2 xl 0 3 ' erg s -1 (~ 10 % 
Eddington for a 1.4 Mg NS) with more than one detector. 
Hence, we cannot determine whether the spectrum evolved 
as the luminosity increased. 

3.4 Period estimation 

Assuming an apparent B magnitude of 28.21 for T14 during 
HI and line-of-sight absorption equivalent to 2.6xl0 21 H 
atom cm~ 2 , Equation 5 yields an absolute V magnitude of 
+ 1.9 for T14 in outburst. The highest 2-10 keV luminosity 
observed from T14 during the outburst was 2.7±0.2xl0 37 
erg s ' 1 for our assumed emission model. We estimated the 
orbital period using Equations 1 and 2. 

We included the uncertainties in Lx, My as well as the 
uncertainties in the relations themselves by simulating 1000 
instances using a Park-Miller random number generator to 
obtain numbers from 0 to 1, and the Box-Muller technique 
to convert these numbers into a normal distribution: if we 
obtain two random numbers a and b between 0 and 1 , then 
y/— 2 In a cos (27r b) and y/—2 In a sin ( 2 - 7 rfe) are two indepen¬ 
dent random numbers with a normal distribution. For each 
parameter, we added a normally distributed random num¬ 
ber multiplied by the ler uncertainty to the best fit value in 
order to obtain the value for that instance. 

For each empirical relation between X-ray to optical flux 
ratio and period (Equation 1 and Equation 2), we obtained 
1000 estimates of the period and ranked these estimates in 
ascending order; we estimated the mean for each period as 
(Pg 40 + L'i 6 o )/2 and the ler uncertainty as (Ps 40 — Pi 6 o)/ 2 , 
where Pi 6 o and Pg 40 are the 160 th and 840 th ranked periods. 
The estimated period from Equation 1 is S.Ol 7 ® hr, while 
Equation 2 yields a period of l.lig '4 hr; these uncertainties 
include the uncertainties in the relations described in Equa- 
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tions 1 and 2 as well as uncertainties in X-ray and optical 
fluxes. 


4 SUMMARY AND CONCLUSIONS 

We conducted a ~13 year programme of Chandra obser¬ 
vations of the central region of M31, monitoring for out¬ 
bursts of X-ray transients. Transients that were partic¬ 
ularly bright or otherwise interesting were followed up 
with HST/ACS, using the F435W filter that approximates 
the B band to search for optical counterparts. CXOM31 
004205.77+411330.43 (T14) is one such transient; the out¬ 
burst was observed by HST in 2011 July, but the follow¬ 
up HST observation occurred only recently, in 2014 June. 
Difference imaging yielded a counterpart with B magnitude 
28.21+0.16. 

Modeling an XMM-Newton spectrum from a previous 
outburst provided a good fit from an absorbed power law 
model with Ah = 2.6±0.3xl0 21 atom cm -2 and T = 
1.84+0.12, with x 2 /dof = 42/54. Fitting the best Chandra 
ACIS spectrum yielded results consistent with the XMM- 
Newton fits. Furthermore, using the XMM-Newton mode for 
ACIS and HRC observations yielded consistent results, sug¬ 
gesting that this is a good description of the spectrum, as 
instrumental differences would result in substantial system¬ 
atic offsets in luminosity if the spectrum was very different. 
Using the best fit column density, and assuming a typically 
flat accretion disk spectrum, we estimated the absolute V 
band magnitude to be ~1.9. 

The highest observed 2-10 keV luminosity during out¬ 
burst was 2.92+0.06X 10 37 erg s _1 with our assumed emis¬ 
sion model. Using two empirical relations between the ratio 
of X-ray to optical flux and the orbital period (Equations 1 
and 2), we estimate the orbital period of T14 to be 3.0j^'2 
hr and l.llg’^ hr respectively. 

We note that Equation 1 was derived from a mixture 
of NS and BH XBs in several different spectral states, while 
Equation 2 was obtained from a more selective sample, in¬ 
cluding only BH XBs in the thermal dominant state. For 
several of our transients, the resulting periods varied by a 
factor >3; however, for others, inc luding T14, the peri ods 
are consistent within uncertainties (IBarnard et alj|2012f ) . 

We also note that relations between optical to X-ray 
flux ratio and orbital period best apply to the peak of the 
outburst, because the optical decay rate is systematically a 
factor ~2 longer than the X-ray decay (IChen et al.11199711 . 
If the peak was unobserved, then the X-ray to optical ratio 
would be higher, resulting in a smaller orbital period. 

It is therefore likely that the orbital period for T14 is 
<4 hr for a BH accretor. This short period may explain the 
rather erratic behaviour during the decay, as shown in Fig.|4j 
because accretion may be modulated by the precession of 
an asymmetric disk. Such behaviour was first observed in 
cataclysmic variables (CVs) with low mass ratios (and short 
orbi tal per i ods). 

lOsakil (Il989h proposed that the super-outbursts ob¬ 
served in CVs with low mass ratios (i.e. short periods) are 
enhanced by a tidal instability that occurs when the outer 
disk crosses the 3:1 resonance with the secondary; the addi¬ 
tional tidal torque causes the disk to elongate and precess, 
and also greatly enhances the loss of angular momentum 


(and therefore the accretion rate). The disk precession is 
prograde in the rest frame, and the secondary repeats its 
motion with respect to the disk on the beat period between 
the orbital period and the precession period, a few percent 
longer than the orbital period. The secondary modulates 
the disk’s viscous dissipation on this period, giving rise to 
the maxima in the optical lightcurve known as superhumps. 
Some short period , persistently br i ght CVs exhibit perma - 
nent superh umps (IPat tersonlll 999l : Iketter fe N avion l2QQ0l '). 
We note that lHaswell et all ( 200ll ~found that X-ray binaries 
with periods <4 hr are likely to have asymmetric, precessing 
disks. 

The X-ray binary associated with the M31 globular 
cluster B158 is observed at a high inclination, and ex¬ 
hibits periodic 10017+50 s (~2.8 hr) variation in the X-ray 
lightcurve in so me observations jTru dolvubov et al.l l2002h 
but not others (IBarnard et al.l 120061 1. The long-term un¬ 
absorbed 0.3-10 keV luminosity lightcurve of XB158 var¬ 
ied between ~4 and ~20xl0 37 erg s , we suggested 
that this modulation was due to c hanges in mass trans ¬ 
fer during the dis k prece ssion cycle (IBarnard et al .1 12013cJ f . 
In IBarnard et al.l (120061 1 we conducted three dimensional 
smooth particle hydrodynamical modeling of XB158, and 
predicted a disk precession period of 29+1 times the orbital 
period, or ~81±3 hr. This motivated us to observe XB158 
30 times over ~30 days with the Swift XRT in order to 
search for modulations in the intensity over timescales of a 
few days. We observed approximately sinusoidal variation of 
the 0.3-10 keV unabsorbed lumi nosity over ^4- 20x 10 3 ' erg 
s- 1 on a 5.65+0.05 day period (IBarnard et al.ll2015lf . 

The outburst of T14 studied with Chandra and HST 
has a decay profile (Fig. [J that is unlike those of Galactic 
transient XBs, which exhibit exponential decay when the 
whole disk is irradiated, and linear decay in case of partial 
irradiation; the decay can change from exponential to linear 
when the X-ray luminosity is no long er sufficient to irradiate 
the whole disk ([King fe Rittenll99a l. We propose that the 
X-ray luminosity may be modulated by variations in mass 
transfer from an asymmetrical, precessing disk. 
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